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Abstract. The size distribution of dust particles in nuclear fusion devices is 

close to the power function. A function of this kind can be the result of brittle 
destruction. From the similarity assumption it follows that the size distribution 
obeys the power law with the exponent between —4 and —1. The model of 
destruction has much in common with the fractal theory. The power exponent 
can be expressed in terms of the fractal dimension. Reasonable assumptions on 
the shape of fragments concretize the power exponent, and vice versa possible 
destruction laws can be inferred on the basis of measured size distributions. 
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1. Introduction 

The dust appears in most of nuclear fusion devices due to the plasma-wall interaction. 
The dust particles pose potential problems to plasma confinement: decrease the 
plasma temperature, absorb tritium, etc. Therefore, it is necessary to investigate 
dust formation mechanisms and parameters of the dust. The size distribution is an 
important characteristic of the dust, as it determines the dust surface area. The latter 
in turn determines the evaporation rate of the dust and the tritium absorption rate. 
To our notion, there is no analytical model explaining experimentally observed dust 
size distributions. 

There are two basic distributions commonly used for approximation of 
experimental results: the log-normal distribution and the power one (the Junge 
distribution) [U [2j [3] . The Junge distribution was observed in several experiments 
O m m [5l [6] for tungsten and carbon dust in the size range from several nanometers 
to tens of microns (Fig.[l]). The exponent (—a) of the power distribution was measured 
to fall between —3.3 and —2.1. The power distribution was observed both in fusion 
devices and in specialized facilities designed for erosion studies. 
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Figure 1. (Color online) Dust size distributions measured at (a) stellarator LHD 
(from ?2]), (b,c) QSPA facility (from (d) tokamak T-10 (from i5j), and (e) 
mirror trap GOL-3 (data taken from \§\). 

The power size distribution was also observed for the atmospheric dust [IJ. The 
distribution was explained by the model of particle coagulation [gj . Coagulation results 
in formation of dust agglomerations which are distributed in size according to the 
power law. In plasma devices, the power size distribution was observed not only 
for agglomerations, but also for single dust particles in the size range from several 
nanometers to several microns |2] . Hence it follows that the coagulation model cannot 
explain the dust size distribution in the whole range of sizes. 
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Four basic mechanisms of dust formation are usually distinguished in plasma 
devices: flaking of redeposited layers, brittle destruction, condensation from the 
supersaturated vapor, and growth from hydrocarbon molecules 1^. If the dust 
grows from the gas (condensation from the supersaturated vapor and growth from 
hydrocarbon molecules), the size distribution differs from the power law [TUj. The 
typical size of the dust formed by flaking is greater than one micron 111] . At smaller 
sizes, flaking can be thought of as brittle destruction. The brittle destruction is thus 
the only mechanism that can be responsible for formation of the observed dust. In 
this paper we show that the brittle destruction can indeed provide the power size 
distribution. 

Brittle destruction was studied theoretically and experimentally [l2llT3l[T4l[T5llT6l 

I17j . The power law for single dust particles was observed in the wide range between 
two typical sizes: the grain size of materials and interatomic distance. This fact 
suggests that the law of material fragmentation is independent of the scale. Moreover, 
the power distribution itself, which intrinsically has no typical size scale, points to the 
same conclusion. In Section [2] we flnd which dust size distributions can be realized 
under the assumption of scale similarity and show that these are the power ones. 
The exponent of the distribution function can vary in the range from —4 to —1. In 
Section[3]we show the interrelation between the discussed model and fractal theory and 
connection of the distribution exponent with the dimension of packing. In Section [4] 
we relate basic features of fragmentation with distribution functions of fragments for 
several representative cases. In Section [5] we analyze available experimental data and 
discuss possible implications of measured power exponents. 

2. Mathematical model 

First we formulate the problem mathematically. The dust particles produced due to 
the brittle destruction are fragments of a solid body. The body is to be divided into 
pieces according to some law. Under the assumption of scale similarity, the division 
law must be independent of the fragment size. Let us mentally remove the pieces from 
the body one by one in the order of decreasing size and number them sequentially by 
the index n. The size distribution function does not depend on the shape of the initial 
body only for fragments of the size much smaller than the typical size of the initial 
body. Hence we assume 1. 

Denote the volume of the body remained after removal of the n-th piece by Vn , 
and its surface area by Sn- We can write recurrent expressions relating Vn, Sn, and 
the characteristic size r„ of the n-th fragment: 



where ci and C2 are coefflcients that depend on the shape of removed fragments. In 
what follows we consider these coefflcients to be independent of the fragment size 
and its number due to the scale similarity assumption. The coefflcient ci is positive 
because the body volume decreases as we remove fragments. The coefflcient C2 can be 
either positive or negative. 

It is necessary to express r^+i in terms of Vn and Sn to close the recurrent 
scheme. After removal of many pieces, the size of next fragments is much smaller 
than typical sizes of the initial body. The structure of the residual body does not 



Vn = K-1 - cirl 

Sn ~ Sn-1 + C2fn 



(1) 

(2) 
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depend any more on the shape of the initial body and is determined by the shape of 
fragments. As n increases, the local structure of the body remains self-similar because 
the construction law is the same at all scales. The difference is only in the size of 
structure elements. From the similarity assumption it follows that the size of removed 
fragments is proportional to the size of these elements. The linear size of the structure 
element can be determined as the ratio of its volume to the surface area. Under this 
definition, the result does not depend on the quantity of elements considered, since 
both the volume and the surface area are proportional to this quantity. Thus we 
obtain the required expression: 

fn+l = C^Vn/Sn, (3) 

where C3 is a positive coefficient dependent on the shape of removed fragments. 

It is convenient to analyze the sequence <?„ = SfjV,^ instead of the sequences 
Vn, Sm and r„. Its recurrent expression found from ([I]), ^ and ([s]) does not contain 
other variables: 

(l + C2ci/g„)^ 

9"+! = In J— 2 • (4) 

This sequence tends to a nonzero constant Qoo or to infinity. The second case is realized 
only under the condition 

3c2 + 2ciC3 > 0. (5) 

In the first case, recurrent expressions (jlj and ([2| can be simplified for n ^ 1: 

Vn+i = V„ (1 - c^cl/q^) , (6) 

Sn+1 = Sn (1 + C2cl/qoo) ■ (7) 

Formula Q in this limit reads as 

{l-Clcl/qooY = {l + C2ci/q^)\ (8) 

It follows from formulae ^ and ([t]) that the sequences Vn and Sn are geometrical 
progressions: 

K oc (l-cic^/goo)", (9) 

Sn CC {l + C2cl/qoo)"' ■ (10) 

Substituting ^ and ^ into ^ gives 

r„ cx (l - ciCg/goo)"^^ , (11) 

where we used expression ([s]) for simplification. Hence it follows the size distribution 
function of the fragments: 

dn 
dr 



fir) 



(x-. (12) 
r 



This result has not been observed in experiments, therefore we focus our attention on 
the second case. 

If the sequence g„ tends to infinity, it is possible to simplify the expression Q in 
the limit of large qn : 

1 + 3^2i+2cic|A ^ ^ (3^^^2 ^ 2cic^) . (13) 
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For n 3> 1, the sequence g„ grows hnearly: 

9n ^ (3C2C3 + 2C1C3) 71. 

The dependence of Vn on n foUows from and ([W]): 

whence we obtain the asymptotic behavior of Vn' 

_ 



(14) 
(15) 

(16) 



Similarly, the expression for the surface area of the remained body follows from 
formulae ^ and ( [HI ): 



The substitution of (16) and (17 1 into (|3| gives 



_ C2<:3+ciC3 



This dependence yields the size distribution function 



fir) 



dn 
dr 



a 



(17) 

(18) 
(19) 



1 + C1C3/C2 ■ 

From the condition ([5|, we find that a falls into the interval from 1 to 4. All known 
experimental results are in this interval. 

In the case of a iV-dimensional body, similar calculations give the allowable range 
for the exponent between —N — 1 and —1. 



3. Analogy with fractals 

The above model of body decomposition resembles an algorithm of fractal construction 
and can be treated in terms of the fractal theory. There is a relation between the 
exponent of fragments size distribution and the fractal dimension of packing. The 
latter is determined as Hausdorff dimension of the residual set [18]. 

The above model of fragmentation assumes division of the body into an infinite 
number of domains filling the body volume completely. The same situation is 
realized at tiling packings, for instance, at osculatory packings by spheres [THl [2U] . 
However, tiling packings are generally non-self-similar. This fact makes the involved 
mathematics significantly more complex. Therefore we derive the key relation with a 
self-similar fractal like Sierpinski carpet [TSl HI] which is not a classical object of the 
packings theory. 

For fractals of this kind, a number of equal-size fragments is removed at every 
construction stage. The process is characterized by two parameters: the ratio of linear 
sizes of fragments removed at consequent stages (fci > 1), and the ratio of numbers of 
fragments removed at consequent stages (fc2 > 1)- Denote the number of stage by i, 
then we can write for the fragment number n and the fragment size r 

n = l + k2 + kl + ... + kl^ ^ — , (20) 

fc2 - 1 



(21) 
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and find the relation between them for i 3> 1: 

n(xr-i°s.ifc2_ (22) 

The logj.^ /c2 equals to fractal dimension D [18j [21]. Thus the fragment size 
distribution function is 

ocr-i-^, a = l + D. (23) 



dr 



The same expression was proved for some non-self-similar fractals, particularly for 
osculatory packings [22l [23]. Possibly it is a general law for every decomposition 
described in these terms. 

The fractal dimension of the residual set in three-dimensional space evidently 
falls between and 3. The interval for possible values of a thus naturally follows 
from the dimension of space. However, the result of Section [2] is more general, since 
fractional-dimensional nature is not assumed there and the exponent is related to the 
shape of fragments. 

Let us illustrate the geometry of fragmentation for different values of a. The case 
of large a may be presented by Sierpinski cube (Fig. [2]), i.e. the three-dimensional 
analogue of the Sierpinski carpet [TH], for which a ~ 3.97. The special feature of the 
Sierpinski cube is that it is "spongy" . 



More rigorously, the combination of coefficients in the expression ( 19 ) can be 
interpreted in this way: 

where AVn — Vn+i — Vn and A5„ — Sn+i — Sn- The difference AVn always equals the 
volume of the n-th fragment V*. The quantity ASn can differ from the area of the 
n-th fragment 5* due to partial coincidence of the fragment surface with the surface of 
the fragmented body. In the Sierpinski cube, the fragments do not touch each other, 
so that ASn = Sn, and C2 has the maximum (positive) value possible for fragments of 
this shape and size. Consequently, the expression ( p4| takes the form 

C2 J* / Sn 

The surface area of the removed fragment is the same as the surface contained in one 
cell of the cube, but the volume of the residual in the cell is 26 times bigger. So the 



value of (25) is small (« 1/26), and a in (19 1 is close to 4 



The three-dimensional analogue of Sierpinski triangle [TS] is a tetrahedron from 
which the central octahedron is removed so that four tetrahedrons of the twice shorter 



edge remain. The exponent in distribution (23) for it is —3 sharp. This is the 
separating case corresponding to C2 = 0, which means the surface area does not 
change in the fractal construction process. If the removed fragment is larger and four 
remaining tetrahedrons are smaller (Fig. [3|, then fci > 2, ^2 = 4, C2 < 0, and thus the 
exponent is smaller in absolute value. The limit of very small remaining tetrahedrons 
corresponds to a = 1, which is another limiting value. 

Note that the discussed fractals do not fully conform to the introduced model 
of destruction: the coefficient C3 oscillates as a function of the fragment number. 
Therefore, it is necessary to take an average value of C3. 

In Ref . [5] , a fractal structure of separate dust particles falling on a substrate was 
detected. The fractal dimension of particles was measured to be 2.2 ± 0.2, which is 
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Figure 2. (Color online) Three generations of fragments (marked by different 
colors) of Sierpinski cube having a segment cut away. The inset shows which 
segment is cut away. 




Figure 3. (Color online) The scheme of fragmentation for the modified Sierpinski 
tetrahedron with the fractal dimension 1.4. 



close to the coefncient a ~ 2.3 ± 0.1 for the same dust. This coincidence conflicts 



with relation (23) that follows from the very basics of the fractal theory. Thus the 



used method of measuring the fractal dimension is open to question, although the 
observation of a fractal structure is important by itself. 
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4. Additional assumptions 

Experiments show significant diversity in values of a which, apparently, is caused 
by diverse fragmentation laws in different experimental setups. It is difficult to 
uniquely derive the fragmentation law from the exponent. Therefore we analyze 
several reasonable fragmentation mechanisms, calculate the parameter a for them, 
and correlate basic features of fragmentation with observable distribution functions of 
fragments. 

First, we assume shape regularity of fragments. In the general case, AVn = V* 
and ASn < S* ■ Whence it follows 



ASn 



> 



V* 

q* 



Substitution of ( 26 1 into ( 25 ) gives 



C1C3 


> 


V* 


1 ^ 


C2 


/ 





(26) 



(27) 



If the shape of each fragment is more regular than the shape of the residual, then the 



right-hand side of inequality (27) is much greater than unity, as is the left-hand side. 
Consequently, the exponent in the size distribution function ( 19 1 differs from —3 by 
a small parameter. In the case of A'^-dimensional body, similar calculations give the 
exponent close to —N. 

Second, we take into account the experimentally observed shape of fragments. In 
Ref . [2] , single small dust particles were photographed, and their shapes are spherical. 
Shapes of large (> 1 \iva) particles are irregular. Filling of space with spheres was 
studied in connection with osculatory packing problem |Tni [S^ . The key feature of 
the osculatory packing is that the sphere of maximum possible radius is inscribed into 
the residual at each step. The osculatory packing in three-dimensional space gives 
a « 3.47 1241 . The osculatory packing by disks in two-dimensional space corresponds 
to a « 2.3. 

Third, we use energy considerations. The brittle destruction is more energetically 
efficient than evaporation, but it also needs energy for surface tension. If this energy 
is a factor, then fragmentation must happen with minimal surface formation. Any 
packing by spheres is energy ineffective since spheres can touch neighbours only in one 
point, and no surface is shared between fragments. An example of energy effective 
fragmentation is the modified Sierpinski tetrahedron for which the surface area of the 
residual decreases and the total surface area of all fragments is limited. 

There is a compromise solution that meets both sphericity and energy 
considerations. At the modified Sierpinski tetrahedron, the parameter k\ is free. We 
can vary it to make fragments as spherical as possible. We take the ratio / as a 
criterion of sphericity; this combination is maximal for a ball. The maximum for the 
modified Sierpinski tetrahedron corresponds to a « 2.4 (Fig.jsj). 



5. Discussion 

Several references [31 [5] give experimental results for a in the range from 2.2 to 2.3 with 
an accuracy of 0.1. Close values are given by osculatory packing by disks (a « 2.3) and 
modified Sierpinski tetrahedron [a « 2.4). This suggests that material destruction at 
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described facilities is either energy-abundant two-dimensional or energy-scarce three- 
dimensional. The former variant seems more realistic since the penetration depth for 
particle energies involved is shorter than the smallest size of observed dust particles. 

There are experiments showing greater values of the coefficient a [Sl[7]. In Ref. [3], 
not only size distribution was measured in the interval from 2 /xm to 40 /xm, but also 
the depth of erosion ~ 10 /im in a single shot. The strong erosion appears due to hot 
electrons that deeply penetrate into the target and explode it volumetrically. These 
results could demonstrate the transition from the two-dimensional fragmentation to 
the three-dimensional one. Indeed, for small particles a « 3.3, which is clearly 
distinguished from the value of 2.3 in [31[S]. Unfortunately, for particles larger than 
10 /xm the statistics is poor and no quantitative analysis is possible. 

The exponent close to —3.47 corresponding to three-dimensional osculatory 
packing by spheres is observed for several non- fusion objects [2S1[5B]. First, the size 
distribution of interstellar grains has the power in the interval from —3.6 to —3.3 for 
variety of materials. Second, brittle destructed materials (coal mine dust, crushed lead 
glass) has the fractal dimension corresponding to a « 3.5. Of course, the measured 
dimension is the dimension of fragment surface rather than that of the residual, but 
they are equal in this case. 

The observed distributions can be modified by collisions of fragments. If collisions 
round corners of fragments like sea waves polish pebbles through their contact, then 
the resulting law of fragmentation will approach the osculatory packing by spheres. 
Thus, the exponent close to —3.5 may bear witness to long-term collisional evolution 
of fragments. 

There are several factors not included into the model which can result in difference 
between theoretical predictions and measurements. Counting the dust particles and 
measuring their size are performed after impaction against the substrate. The particles 
can coagulate or change shape on the way to the substrate or on the impact, can 
overlap each other on the substrate, etc. Therefore, to extract maximum information 
from the developed model, a great care is to be taken to minimize evolution of the 
dust on the way to the detector. 
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